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Abstract 


Many new Earth remote-sensing instruments arc embracing both the advantages and added com- 
plexity that result from interferometric or fully polarimetric operation. To increase instiumeut, 
understanding and functionality, a model of the signals these instruments measure is presented. A 
stochastic model is used as it recognizes the non-del orininistie nature of any real-world measurements, 
while also providing a tractable mathematical framework. A stationary. Gaussian-distributed model 
structure is proposed. Temporal and spectral correlation measures provide a statistical description 
of the physical properties of coherence and polarization-state. From this relationship, the model is 
mathematically defined. The model is shown to be unique for any sot of physical parameters. A 
method of realizing the model (necessary tor applications such as synthetic ralib] at ion-signal gen- 
eration) is given, and computer simulation results are presented. The signals are mist meted using 
the output of a multi-input, multi-output linear filter system, driven with white noise. 


1 Introduction 


Radio-interferometers and polariinoters arc being list'd mure widely in remote sensing toi probing tlie 
Earth’s lands and oceans (e.g.. (Iven. 1998. Piepmeier and Casiewsld. 2001].) While these techniques 
have a rich histoiy in the space sciences, Earth-viewing systems an' relatively new. There are 
distinct, differences in the operation and calibration techniques of space-viewing telescopes versus 
orbiting Earth-viewing instruments. Operationally. Earth-viewing imagers haw a fraction of the 
integration time available compared to radio-telescopes: seconds vs. hours. Additionally, radio- 
telescope elements have extremely narrow beams (~ 0.5'’) compared to low-Earth orbiting imagers, 
whose wide beams (~ 00 90") are required to cover a reasonable swath. Perhaps more important 

are the calibration differences. Interferometric radio-telescopes are typically calibrated using extra- 
solar point sources, whereby the system can he characterized to within a common gain coefficient 
[Thompson et al.. 1991]. Furthermore, if the flux of the point source is known, the interferometer ran 
be absolutely calibrated. There are no obvious point sources when looking downward at the Earth, 
at least ones that emit energy within the protected spectrum for passive observing. Fuithennoie. 
because of their large beaimvidtlis, Earth- viewing interferometers cannot selectively view a single 
extra-solar point source. (Besides, the time and operations required for regularly rotating a spacecraft 
for calibration would cause data loss and increase mission costs.) Therefore, a different calibration 
technique must be devised for orbiting Earth-imaging systems. This paper lays the groundwork for 
such a technique by rigorously examining the signals that these radiometers measure. 
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While polarimetry and interferometry arc usually investigated independently, they are based 
upon the same idea measuring the interdependence of two signals. For polarimetry these signals 
are the vertical and horizontal field amplitudes (or some equivalent pair.) Two-beam interferometry 
involves measuring the coherency of two signals separated in space and/or time. In older to fully 
understand these instruments, it is desirable to have an mrnr.it e model of the types of signal pairs 
they measure. This type of model can lie employed in mathematical analysis, computer simulations 
and for the generation of synthetic calibration signals (one possible on-orbit calibration tool.) 

Electromagnetic, waves generally have some degree of randomness an nnpolarizod component 
in polarimetry or an incoherent component in interferometry. For this reason it. is logical to employ 
a stochastic method when modeling the pair of signals. The relationship between the statistics of the 
model and the standard terminology in polarimetry and interferometry should be well understood. 

The model structure, statistical parameters and constraints art' described in Section 2. Sec- 
tions :i and 4 describe how this model s temporal and spectral correlation hmctions can define 
poiarimetric and interferometric properties. Using this framework. Section r, shows how a com- 
prehensive model can be uniquely determined from the physical propcrtii's of the modeled wave. 
Section ti gives a computational method for realizing the model. A Summary is then presented 
followed by an Appendix containing mathematical derivations. A polariniotiic example is carried 
through Sections 3, 5 and (> in order to demonstrate the ideas presented. 


2 Model Structure and Basic Properties 


2.1 Mathematical Signal Representation 


The standard literature in both polarimetry and interferometry makes extensive use of a 'narrowband* 
signal representation (e.g. [Hcrilt, 1987. Brosse.au. 1998].) This representation alludes to the fart 
that the measured signals exist within a finite frecpieuey hand (eentered on some frequency x*n.) The 


genert 


>ral form of this representation is shown below. 


/(/) — f/(t)eosu.'ot /j(t)siiw () t 


= He{ [«(/) + ib[t) 

It should be noted that a narrowband representation does not imply a small fractional bandwidth. 

As mentioned in the previous section, the model must describe a pair of signals. These signals 
will be labeled X(t) and Y{1) and have a narrowband definition as given below. 


X(t) = ,l(t)msii)t B(t)siiu'of 

= R e(P(/ 

( 2 ) 

Y(t) - <*'(f)eosu.\if f?(f)sin^iif 
- Re{Q(f 

The functions are capitalized to indicate that they are random processes. Il will also be assumed that 
the model is stationary and Gaussian distributed. Stationarity is a common assumption in optics 
[Ilecht. 19871 and a Gaussian amplitude distribution can be justified by appealing to the Central 


Limit Theorem. The consequences of not assuming a Gaussian distribution are discussed at the end 



of Chapter 5. 
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In order to ensure stationarity in the above signal model (Equation 2), X{f) mid 1 (t) must 
bo jointly-st.atkai.ir>-. For simplicity it. will bo assumed that -1(0- B(t). C(t) and D(t) (and hence 
P(t) and Q(t)) must, also be joiutly-stationary. To produce a Gaussian-dist.ributed model for X(t) 
and Y(t) it will be assumed that -1(0- B(t). C(t) and D(t) are Gaussian processes, which leads to 
l)ivariate Gaussian distributions in the complex plane tor Pit) and Q(t). 

The random processes A{t). B[t). C(t) and D(t) air real, joint ly-stat ionary and Gaussian. This 
means their multivariate probability density function is completely defined bv their means and second 
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order moments (shown below.) 


E[A(t)\ - 


E[D{t)} = hb 
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Rr(r) 

E\D(t)D{t 
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E[A(t)B(t 

D] 

~ Rad(t) 

E[A(t)C(t 

t)\ 

-- R \c{~* 

E\A(t)D(t 

t)S 

- R-ad( t ) 

E\D{t)C(i 

r)] 

-- Rdc'(t) 

E[D{t)D{1 

T)\ 

- Rud(t) 

E\C{t)D{t 

r)] 

- Ecd(t) 


Th<> functions above can be used to create a covariance matrix and mean vector for any set of points 
within the signal. As these two properties define a multivariate. Gaussian distribution, the above 
functions give a complete description of the model. 


It is also possible to find the means and second order moments for P{1) and Q(t) b\ noting that 
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P(t) = A(t) 4 i B(t) and Q(t) = C(t.) 4 i D(t.) (see E<iuatinns 1 and 2.) 

t , P - E\Plt)]=r-E[A(t)AiDlt)\ = f i. x + inB 
,i Q - E[Q(t)} - E\CV) + i D(t)} = nc 4 > I'd 
Rp{t) - E[(A(t) 4 i B(t))(A(t - t ) i B(t - r))] 

— R a (t) 4 R[i(t) - HR.\b(t) - R,\u( T )) 

(4) 

Rq(t) - E\(C(t) 4 i D(t))(C(t t) i D{t r)'l| 

Rr{~) 4 Rd(t ) - HRcd(t) Real- t)) 

Ri>q{t) - E\(A(t) - iB(l))(C(t r) i D{t t))\ 

— R U'(") 4 R[ia( T ) i(Rad( t ) Rbc{t )) 

2.2 Correlation and Spectral Functions 

It. is often useful to work with these; second order moment, functions in the Fourier domain. Tlu*se 
frequency domain representations are known as spectral functions and will he denoted b\ an .S 
i.e. A{Rvv(t)} - S rv {u.’). where d{ } denotes the Fourier transform. 

Correlation functions and their corresponding spectral functions cannot take on any arbitrary 
form there are theoretical restrictions on the range of legal function sets. These restrictions 
(given below) are most easily presented with respect to the. spectral functions. 

87 du-’) > 0 .SV(a>) > 1) .SV r(-<j)| £ v' (5) 

It can be shown that a complex, bivariate, stationary random process can be generated to have any 


spectral functions that satisfy Equation - r >. 


y 


It was stated in the previous section that, a real, stationaiy. multivariate. Gaussian ])roeess is 
uniquely described by its means and correlation functions (or equivalently spectra.) This is not the 
case for a complex, stationaiy. multivariate. Gaussian process. This can be demonstrated by exam- 
ining Equation 1 and noticing that there are multiple sets of real and imaginary components (-'1(f). 
D(t). C(t). D(t)) that will result in the same correlation functions tor P(t) and Q(t) i.e. speci- 
fying the correlation functions of P{t) and Q(t) does not uniquely define the model. To completely 
determine the model it, is necessary and sufficient to define the functions given in Equation If. 


3 Polarimetry — Background and Relation to the Model 

This section presents published material [Paneliaratnain. 1075. Brosseau. 1998. fshimaru. 1991] in a 


way that relates it to the proposed model structure. 


3.1 The Stokes Parameters 


Polarimetry involves studying the relationship between the horizontal and vertical fields of an elec- 


tromagnetic wave. These two fields arc usually written in narrowband fonn. 


Eji(t) --- Re{E,(f 

E v {t) - R ('{^(fjc 1 -') 


( 0 ) 
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The polarization state is generally measured using the ; Stokes parameters’. These parameters are 
defined in relation to the narrowband formulae above. 


- 


- 

1 


< |£i(0| 2 > + < |E 2 (t)| 2 > 

Q 


< |£,(/)| 2 > - < |£ 2 (0| 2 > 

V 


2 !{('.{< E l { 1 )Er,(t) >} 

V 


2Iin{< Ei (1)Ei {t) >} 


This indicates that four parameters are needed to uniquely specify the polarization state of a wave. 
These four parameters are known to obey the inequality given below. 

I > y/Q 2 + U- + V* (*0 

The ratio of the two terms in Equation 8 is known as the' ‘degree of polarization . 

Earth sensing radiometric instruments lend themselves to the use of the ‘modified Stokes parame- 
ters'. These are typically used when the vertical and horizontal polarizations correspond to preferred 
axes tangent and normal to the Earth s surface. 


*1 


< |E!(f)| 2 > 


itfl 

2 

s> 


< I E 2 {t)\ 2 > 


LJ1 

2 

s-> 


2Rp{< Ei(t)£i(0 >} 


u 

*1 


2Iin{< £i(0E.?(f) >} 


V 


The parameters *1 and -s 2 represent the power in the horizontal and vertical channels respectively. 
In the past many Earth-sensing instruments have only measured these two properties. However. 
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remit- Earth-sensing instruments art* now measuring the full set of four parameters this in< teases 
the need for a general model of polarized signals. 

The inequality in Equation 8 ran easily be rewritten in terms of the modified Stokes parameters. 


The result is given below. 



(10) 


It is also obvious from the definition in Equation !) that the inequalities below are true. 


*1 > 0 *2 > 0 
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The polarimetrie theory given above can be related to the model struetiue by invoking stntion- 
arity. Statiouarity implies that the. probability of any given polarization state is constant across the 
signal. In this ease, statiouarity also implies crgodicity which means that the expected values of the 
random processes are equal to the corresponding time averages. If X(f) is set to model Ejj(i) and 
Y(t) models E'v(t) then the modifies! Stokes parameters can be related to the correlation functions 
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as shown below. 
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2R«{ /?/>(?(»)! 

2Iiu{i?rg(0)} 

j 

./ Sq(x)<Iu! 

2R<‘{ J SpQ(u.')d-j} 

2lm{ j S P qU)<E’} 

3.2 Frequency Distributed Stokes Parameters 

The analysis above shows that the polarization state of the modeled wave is determined by the 
area under the spectral functions. While this is correct, a more general interpretation is given m 
Panrharatnam [1!)75]. In that paper it is shown that each frequency component of the wave can 
be treated as a single, independent wave with correlation properties dictated by the- spectra at that 
frequency. The correlation functions an; a measure of the interdependence between or within the 
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signals. As the spectral functions show how this interdependence is distributed across frequency. the 


polarization state' at each frequency is determined by the spectral functions 


- 


" 



.S'p(u. ! ) 

• s 2(u0 


Sq(uJ) 



2R e { Sp(j (u.- ) [ 



2Im{.S> Q (u;)j 


It should be uoteel that the variable v above re'prese'iits a measure of distance trom the central 
fmptency ^ tl . Also, the Stokes parameters in Equation 13 are actually densities i.e. their units 
tire ])ow(T per unit bandwidth ratlin* than just power . 


To get the final wave, the components are considered to be incoherently (as t hey mr. independent) 
summed over frequency. This is intuitively satisfying as spectral functions sum under independent 

addition and Equation 13 becomes Equation 12. 

If Equation 13 is substituted into the inequalities given in Equations 10 and 1 1 tin'll the following 

inequalities are produced. 


2 \/ Sr (i*j)Sq{u>) > 2 \f R e { Spq (u; ) } 2 4 Im { S pq ( ) } ‘ 

£ \/ Sp(uj)Sq(u>) 


Sp(jj) > 0 Sq(u!) >o 

It can be set'll that these art' the same inequalities as those stated in Equation 5. This means that 

the' theoretical restrictions on the. model given in Equation 5 are equivalent- to restrictions on the 
modified Stokes parameters and lienee do not limit the range of polarization state's that can b< 
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modeled. 

The conclusion that can ho drawn from this section is that, the polarization state (distributed over 
frequeney) to be modeled can bo used to define the model parameters S p (jj). Sq(ui) Spq(ui) 
(which in turn can be used to define B p (t). Bq(t) and R P q(t).) This is demonstrated on an 
example signal with a contrived, frequency-dependent polarization state. 

Example: The signal to hr: modeled lias 5 distinct polarization hands within its overall bandwidth 
■>B\V . The lowest fifth- of its spectrum is tmpolanzrd and has unity power density ( modified Stokes 
densities' are [0.5. 0.5,0. 0j^; the second fifth has a degree of polarization of 0.5. unity power density 
and the polarized component is left-circular ([0.5. 0.5. It. 0.5] j; the central fifth is completely linearly 
polarized at 30" and has unity power density fl.75, .25, 0.80(5, 0]A‘ the fomtli section has a degne a .f po- 
larization of 0 . 5 . unity power density and the polarized component, is right.- circular ([0.5. 0.5. 0, 0.5] J: 
and the highest fnqueneg section is also unpohirized with unity power density ([0.5. 0.5. 0. 0 ] ). By 
using Equation 11 the model parameters shown in Table 1 tm found. 

The modified Stokes pammeters for the total wave can he found by performing the integmtums in 


Italics used 
for example 
text 


Equation 12, This rrsults in the modified Stokes vector . 0.<)Z?lk O.MhBH'.Oj. 
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4 Interferometry — Background and Relation to the Model 


4.1 Basic Coherence and Interference 

Two beam interferometry measure* interference and is closely related to the idea of coherence, as 
explained in [ITeeht, 1087]. Interference is produced by tin* addition of two waves that have been 
separated in time and/or space. It is generally assumed that polarization effects can be ignored 
and that the space/ time separation produces a time delay t between the signals. The time delay is 
typically produced by a difference in propagation length ami hence in time of flight. Stationarity is 
also usually assumed, so the two signals cun he modeled as A (/) and 1 (/ — r) as given in Equation 2. 
Interference occurs when the two beams art! addl'd. This sum is expressed below. 


Z(i) - X(t) + Y(t - r) - R o{[P{t) + Q(f - T)c"^ r \t^ nt ) U r >) 


It is the time- averaged, intensity of this sum (the square magnitude of S(t)) that is observed in 
interferometry. Constructive interference occurs when P(t) and Q(t — r) art' in-phase and add to 
produce a large magnitude. Destructive interference occurs when they are out-of-phase and add to 
produce a small magnitude. This gives bright and dark interference fringes respectively. 

/ - <Z{t)Z'{t)> 

- < P(t)r* (t) > I < Q(t r)Q*{t r ) > 


+2Ro{<P{t)Q"(t-T)r u ' nr >} 


(Ki) 


rn(0) + r.j.,(()) + 2Re{r 12 (T)} 



The final step makes use of ergodieity. The filia tions Ta(r). r 22 (r) ami Tiafr) are known as the 
self- and nmtmil-eoherenee. functions. They are eoinnamly used to define the eoherenee of a pair of 
waves and are defined below. 


r u (r) - 

< P(t)P*(i 

t) > c i * ,,T - Rp(t)> ‘~" t 


r, 2 (r) - 

< Q(t)Q*(t 

7) > - Rq{t)c‘^‘ t 

(17 

r 12 (r) = 

< P(i)Q‘(t 

t) > = Rpq(t) r^ r 



Equation 17 shows how the coherence properties of Mu? signals to be modeled determine the model 
parameters R P {t), Rq(t) and Rpq(t) (which in turn can be used to define S P {*), Sq(sj) and 

$pq(^)') 


4.2 Useful Functions in Coherence and Interference 


There are a variety of useful measures defined by the coherence Inactions a brief summary of 
some of them is given here. 

The final term in Equation 10 is known as the interference term. Mathematically, it is this term 
that produces the interference fringes (the other two terms are simply constants.) So. the shape 
of the fringes governs the form of Ti 2 (t), which in turn determines the model parameter Rpq(t) 
(and Spq{v).) The model parameters Sp&) and Sq(jj) are the power spectral densities of the two 
signals i.e. they determine how the power of the signals is distributed about Also, if a signal 
were to interfere with itself, it would he the self- coherence functions (given by R.p(r) or Rq{t)) that 


would produce the interference. term. 
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Ail important interferometric property is the fringe 'visibility function . This function is defined 


bv the intensity given in Equation 15. 

I, nor In, in _ 2|f? rq (-)l 

“ i, no,- + " r„(0) + r 22 (0) “ R P ( o) t RqO>) 


(18) 


Here it is assumed that the oscillatory exponential term in Ti 2 (t) varies much taster than the func- 
tion's magnitude. It is often the case that r u ( 0 ) - r 22 ( 0 ) and the visibility function in Equation 18 
reduces to the absolute value of the •complex degree of coherence' (V I'Via (t) |). 

, r 12 (r) R poiry^ mn 

' ' v /T77IiT)T' 2 ^(f) v/'Mipi 

A similar function to the complex degree of coherence is the ‘complex degree of spectral coherence 
(see [Goodman, 1(185].) This function gives a measure of the distribution of the degree of coherency 
across frequency (a similar idea to the distributed polarization concept presented in Section 3.) The 
function is expressed in terms of the spectral representations of the coherence functions (f?n(u/) — 


jfrnW), G- 22 M = o{r,,(T)} and Q n {*) = d{r 12 (-)}.) 


l> 12 (-t) ~ 


Q 12(^) 


SpojoJ ^-'n) 


1 1 Mfe (wC ) y & / ’ ( ‘-o' — ~ a-'o) 


(20) 


The following properties of the spectra can be shown to be true. 


£uM>0 fc 2 (u?)>n o < |// 12 u-)j < 1 


( 21 ) 


The above equations imply that the inequalities in Equation 5 are always satisfied. This means that 
the theoretical limits on the model do not limit the range of partially coherent signals that can be 


modeled. 
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This section has shown that there' an' a variety ot related interferometric and coherency properties 
(e.g. tlie coherence functions or the interference term, the complex decree ot coherency. the visibility 
function, etc.) that can be used to define 1 the model. As in the polarimetric case, the functions Sp(u;). 
Sq(jj) and Sr Q (uj) (or Rp{r). Rq(t) and Rpq(t)) can all be specified by the physical property's of 
the wave. 

5 Completely Defining the Model 

To completely define the model the functions in Equation 3 must be specified. In Sections 3 and 1 
it was shown how the physical properties of the wave's to be modeled (in either the* polarimetric 
or interferometric cased can be used to define the' correlation functions of P{t) and Q(f). Ilowcvei . 
in Section 2 it. was shown that these correlation functions are insufficient tc> completely define the' 
model (as there are many sets of the functions in Filiation -1 that will result in the 1 same correlation 
functions for P(t) and Q(t).) 

It can also be shown that while exte*nsive use' has been made of the 1 fact that P(t) and Q{i) 
are stationary (as are /1(f). B{t ). C(f) and D(t)). this does not guarantee that A(f) and V (f) are 
also stationary. In this section (and the corresponding Appendix) it will be' shown how enforcing 
stationarity on X(t) and Y(t.) allows the model to be completely defined using only the' correlation 


fiuietions of P{t) and Q(f). 
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5.1 Consequences of Enforcing Stationary 

The.' first step in enforcing stationarity is to give the constraints it places on the* functions in Equa- 


t.ion 3 (see Appendix A.) 


PA — ft B ~ PC ~ I'D 


= 0 


( 22 ) 


Ba(t) 

- Rb(t) 

Rc(t) 

= Rd(t) 

Rac(t) 

— Rbd( t ) 

Rad(t) 

- - Rbc(t) 

Rad{t) 

- Rad( t) 

Rcd(t) 

— R('d{ 


(23) 


The above equations are important because they show that rather than having four constants and 
ten functions to find (as in Equation 3), there are now only six functions (two of which must be 
odd-symmetric.) 

An interesting result can be stated from Equation 2‘h By looking at a single signal (say X(t)) 
and its in-phase and quadrature components (A(t) and D(t)) it can bo. seen that RaW “ RdW au( 1 
R.\b{U) = 0. This means that at any given time the in-phase and quadrature components have the 
same variance and are micorrelated. Since this is a Gaussian process, the result is that the amplitude 
of the signal is Rayleigh distributed and the phase is uniformly distributed. This result is stated m 


[Brosseau. 1998], 
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Equation 23 also results in X(t) and Y(t) bring zero moan and having a well defined set of 
< nit ‘la t io n f i met ions . 

R\ (r) — Ra{ t ) coSv^ot T Rad ( r ) sinu. ? o r 
Ry(r) Rc{~) eosu;oT T Rc;d(t) sinu^r 
/?xv(t") “ Rac(t) T /?ad( ?") siller 

By substituting the relations in Equation -2 and 2d into E(|uat.ion 4 it can he seen that P{t) and 
Q(t) must be zero mean and have correlation functions as given below. 


n P (r) 

= 2 (i?. 4 (r) - 

i Rab(t)) 



- 2 (B c (t) 

i R<:n(~)) 

(25) 

Kpq{t) 

= 2(R M :(t) 

- i Rad{t)) 



The equations above are significant because there is now only one real function defining each of the 
real and imaginary parts of the spectra of P(t) and Q(t) (c.f. Equation 1.) 


5.2 Mathematical Model Definition 

Equation 25 implies that there is now only one set of correlation functions tor A(t), B(t), C(f) and 
D(t) that will produce, a given set of correlation functions for P(t) and Q(t). This set of correlation 
functions is found by simply taking half the real and imaginary parts of Rp(r ). Rq(t) and Rpq{t). 
This results in the six functions required in Equation 2d. It should be noted that Equation 2.Ts 
odd-symmetry condition on R.\n( T ) and Rcd( t ) i s always satisfic’d as Rp(r) and Rq(t) aie always 
c *oujuga te-svmi net ric • (this is a ( onsequenee of Sp(jj) and being real, as im])lied in Equation 5.) 
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This method of determining the correlation functions of .1(f). B(t). C(t) and D(t) c an also be 
expressed in tho frocpioncy domain. 


Ra(t) = 

r t Rp{Rp (7")} - \{Rp(t) + R}>{t)} 

=> Sa(&) ~ 

\ [Spl^j) f Sp ( *- u;) ] 

Rc(t) = 

= \{Rq(t) \ R^(t)} 

s (•{■*') = 

(u.') + S Q (-<*•)] 

Rad (t) 

^Im {J?p(t)J = \{Rp{t) Rp(r)} 

Sad(^) ~ 

j[Sp(u;) Sp{ w)] 

Rcd( t ) — 

^Im{7 ?g(r)} = i{R Q (r) - R^{t)} 

Scd( jj) 

■j[Sq M Sq( ^)] 

Rac(t) - 

±Rc{Rpq(t)} - \{RpQ(r) + Rpq(-)\ 


\1$PQ{^') -t- SpQ(-iO ) j 

Rap( t ) ~ 

-ylni{7?pf v )(r)) - -j{/?PQ(r) 

=► $adM — 

ftSpgM - S^(-w)] 


By using the equations above and the first four equalities in Equation 2d. all the functions in 
Equation It can be found. This leads to a model that is completely specified; and, that specification 
is unique (i.e. it is the only one that will produce the desired correlation functions for P(t.) and Q(t) 
while maintaining stat.ionarity of A(f) and 1 (f) ) 

It can be seen that tin' even- and odd-symmetric portions of the. spectra -VMr imt * 

S PQ (u!) are being taken. This results in the above spectra being con. jugat (-symmetric as would be 


expected for the spectra of real, random processes. If the complete set of spectra are arranged to 
form the spectral matrix, it can be. shown that, the result is a Ilermitian. positive semi-definite matrix 
(see Appendix B.) This is a general property of all spectral matrices (see [Jenkins and Watts. liKiK]) 
and in the two-procoss ease gives the inec qualities of Equation 

Example: This process will hr applied to the spectral [unctions from the example in .Section J. 
Using Epilation 26 tin’ results shown in Table 2 are found. 

5.3 Summary of Model Design 

A brief summary of the model design is presented in the tot in of a list. 

1. Define the spectral or correlation functions of P(t) and Q(t) by considering the physical prop- 
erties (either interferometric or polarimetrie) of the waves to be modeled. 

2. Apply Equation 2(i to get. six valid spectral or correlation functions tor .4(f). B(t). ( (t) and 

D{t). 

3. Apply Equation 23 to get the complete set of statistical functions (as given in Equation 3) of 
the random processes A{t). B(t), C(t) and DU) and thus completely define the model. 

It should be noted that there is no loss of generality in any of the steps, and that a single set of 
physical parameters can only produce one valid model. It can also be seen that if the Gaussian 
assumption is not made, then the model is not uniquely defined (as second order statistics no longei 
completely describe the p.d.f.). However, the previous analysis can still be used to define a stationary. 


covari aiioe-ergoc li<* model u]> to second-order statistics. The Gaussian model is a special ease' of this 
where stationary implies ergodicity and second-order statistic s specify all higher orders. 

6 Realizing the Model 

The previous sections have developed a stochastic model for a pair of signals. This model is defined hv 
the temporal and spectral correlation functions ot A(t ). B(t ). C(f ) and D(t). These ran be found from 
physical signal properties and used to specify the probability density functions of the model. This is 
a clear mathematical model that has the potential to be useful in theoretical analysis. Howcvei. toi 
applications such as computer simulations and synthetic signal generation it is necessary to create 
realizations of these random processes. It is sufficient to create realizations of the real variables 
.4(f). B(t), C(1) and D(t). as the other variables (P(f). Q(i)- A r (f) and \ (f)) ran be created by 
deterministic functions of these four realizations. 

As the probability density functions are known, it is possible to e rente, a sampled realization 
direc tly. If N points of data were required, the correlation matrix for these points could be calculated 
and a multivariate. Gaussian random number generator applied (such as the mvnrnd command found 
in the MATLAB software' package.) However, difficulties arise when the dimensionality of this p.d.f. 
is considered. There are four co-dependent outputs at every sample so a signal of length N would 
have a p.d.f. of dimension 4 N. For example, a signal of 125 MIIz bandwidth requires a minimum 
sampling rate of 250 x 10 r> samples per second thus a two-second signal would require 500 x Hf 


Emphasis 
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software 
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samples and result- in a probability density function ot dimension two billion. This is computationally 
impractical so a simpler method must bo tound. 

A common technique in one dimension is to use a linear filter to shape noise into a desired 
spectra! shape. The problem here is more complicated as the function is from one dimension 
(time) into four dimensions (A(f), B(t). C(t ). D(t)}. A generalized filter stiurtme with A in- 
puts I,v(f)) and M outputs (Oi(t).(h(t) 0^(1)) is given below (this structure 

is presented in [Jenkins and Watts. 1998].) 


Oi(t) = 

V A 

L^.) ■ 


02(0 - 

E; v 

-jMO* WO 


OM{t) =• J2r 1 * ^(0 

The filter responses must he chosen so that the desired correlation functions are realized. These 
functions are given by the expressions below [Jenkins and Watts. 1908]. 


Ro t .o,M) - E[O p (t)O n (t - r)j (- Ro q o v (- r)) 

= E [(E; V 1 M*> * EH)) (Ef , M* r) * Ij(t T))j 

- E [E;-1 E;-i (M'> * E(t)) (h qj (t - T) * /,(/ r))] 


= ET Ef-i (0 * V H) * K ( ~r) 

This set of equations has a simple]- form in the Fourier domain, as the. convolutions become multi- 
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plications. 


EihEiLi 


= E; v i E; v i s,^)H,Ju:)h; i:i (*) 

The final stop takes advantage of conjugate symmetry (which is guanmtoed by the fact that 
is real for all mn.) 

A general sot of solutions for the system of equations given by Equation 29 is lion-trivial as tlio 
syst(*m is non-linear. Simplifications can be made by making certain assumptions the first of 
which is that the input random processes are independent, white. Gaussian, unit variance, and zero 
mean (so S\,x } M = £(/ -./).) This ensures the outputs will be Gaussian and zero mean as required; 


and that Equation 29 reduces to the form given below. 


So„oM) « Ei V iE; j)n tl> U)H^) ( 

- e; v i 

In this ease only four output processes art* required (A/ ™ 4.) This means we have' ten independent 
equations one 4 for each of the correlation functions given in Equation 4 (there are lli equations but 
there is redundancy, as shown by tlie bracketed term in Equation 2K.) A solution is piesented for the 
case of four input processes (i.e. A ~ 4) although the nic*tliod ean be applied to larger dimensions 
provided that M = N and the specified spectra satisfy the properties of a spectral matrix. The 
equations for this problem (as defined by Equation 40) are given below (the us dependence has been 
dropped for clarity.) The outputs are denoted by A. D. C, D as before, while the inputs art' numbered 
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S A - |//.u| a + |//.d 2 + l^.i| 2 + |//.4i| 2 


Sab 

il 

&S 

+ + H-uHh, 

Sac 

= H M HZ: 1 

*f II \2 H C2 h IIaaIIca T 

Sad 

= iiMiib 1 

T II A'2 II D 2 "f HaaIIpa ^ II A [lip 1 

Sb 

= \H m \‘ + 

|tffej|' 2 T |i/r?;d 2 + \Hbi\ 2 

Sue 

- UmlllA 

t II D 2 He "2 "1” II B :\II c i d Hu\He^ 

S BO 

= IbJIbi 

•f II B2 II D2 "I II BA H PA I II Bill pi 

Sc 

= l^cil 2 + |H r2 | 2 + |//o| 2 + l^cil 2 

Sep 

- HciUpi 

\ Ilcilip 1 //r : t /e n f fla Jip 

So 

= |J?£>l| a +|ffM| a + I^OT| 2 + I^Zw| 2 


Although these equations are non-linear they can still he solved relatively easily. The solution method 
developed is most sueeinetly expressed in matrix notation. 

H , u Ha‘2 Hm H A4 


II - 


II B 1 ?! B2 H BA II BA 


( 32 ) 


H C i He > He II ci 


Hp 1 Hpn Hp-,\ Hpi 
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S A Sab Sac Sad 
S*\b Sb Sbc S B d 

(33) 

Sac She Sc Sen 

Sad S bd S(* d So 
Using this notation. Equation 31 can bo rewritten a* shown below. 

f^77 77 + (34) 

For mat rices the 1 opeTation represents a conjugate transpose. In order to find a suitable sot ot filters 
it is sufficient, to solve Equation M at the frequency points of interest It is shown in [Strang. 197(>] 
that because S is ] >osit i ve semi- definite (see Appendix. B and [Jenkins and Watts. 19(i8]). Equation d 1 
can always bo solved. A method for doing this is outlined below. 

Because the' spectral matrix $ is Honnitiaii. it is always diagonalizable, its eigenvalue an' real 
and its eig('iivectors are orthogonal. This allows a simple* solution for Equation 34 to be found In 
using the diagonalized form of 5. If the eigenvalues are arrange'd on the main diagonal in the matrix 
A and the* eigenvectors in the matrix E then the following is true. 

f 

S - !vava! [ 

(35) 

S = (f VaV 

f = (f VA)(I v / A) t 

\/a is a matrix with the square-roots of the eigenvalues on the main diagonal (so that V A \/a — A.) 




Positive* semi-definiteness giuuantees the oigenveetors are positive, so V A will always he a real matrix 


and hence he (‘qua! to its own conjugate transpose. 

By comparing Equation 35 with Equation 34 it can be seen that the filter responses can be 

calculated by as follows. 

77 -- !\/X (36) 

Equation 30 shows how Equation 30 can be solved by finding the eigenvalues and eigenvectors of 
the spectral matrix *5 at each frequency point. It is easy to show that this method will produce 
filter spectra that are conjugate-symmetric, which is necessary to ensure that the filters have a real 
response. Reordering the eigen- values and, -vectors will still result in a valid solution at a single 
frequency point but care should be taken not to do this when constructing functions over many 
frequency points. Doing so would result in a sharp discontinuity in the filter spectia pi od need wliic h 
would increase the length of the filter impulse response. 

In this section it, has been shown how independent, white. Gaussian noise functions (which aie 
easily generated) can be fed into a system of linear filters to produce realizations of the model. The 
filter sets are generated using Equation :J(i which depends on the model parameters. This process is 
computationally tractable and produces results like those shown in the example below. 

Example: The example functions used in .Sections 3 and 5 were realized using 1 the methodology* 
given in this section. This process was carried out using MATLAD and plots of the resulting spectra 
can be seen in Figure L It can be seen that the frequency axes extend outside the region [ B\\\ B\V] 
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this corresponds to ovevsam j ding. The free/ i h 'ncy axi \s an > 1 iorn in livx d si ) that th e range ot uunlinsed 
fcaiuencios fulls between —0.5 and 0.5. 

The estimates of the realized spectra were found using periodogram averaging. 100 spectra were 
averaged in each case arid a rectangular window was used in the time domain to remove the noisy 
terms associated with a large |r|. Each plot, has M) frequency points and the filters were truncated 
to 199 taps . Figure 1 shows that the resulting spectra agree closely with those specified The small 
differences can be accounted for by the necessary truncation of the generation filters and by the fact 
that a finite number of spectra were list'd to create tin 1 periodogram average. 

7 Summary 

This paper begins by showing how standard measures in polarimotry and interferometry (such as 
Stokes parameters and coherence functions) can be interpreted in terms of the model s temporal 
or spectral correlation functions. This statistical interpretation then allows a stationary, Gaussian 
model of the signal pair to be defined. A stationary. Gaussian model can be justified physically. 

This model structure was shown not to limit the range 1 of physical signals that can be modeled (i.e. 
all coherency- or polarization-states could be modeled.) Additionally, the stationarity assumption 
was shown to lead to a unique model for any given set of physical properties. This indicates a 
comprehensive, well defined stochastic model in either the polarimetric or interferometric paradigm. 

To realize this signal (as would be required in such applications as synthetic* calibration signal 
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generation), it is necessary to produce a signal pair that has the same properties as the model. This 
can be done by generalizing a well-known noise-shaping technique in which a white, Gaussian process 
is passed through a linear filter in order to color its spectrum to a dt'sired shape, The method was 


corroborated by presenting results from a MATLAB implementation. 


Appendix A Derivation of Equation 22, 23 and 24 

Because X(t) and Y(t) are Gaussian, they will be stationary if mid only if their means and second 


order moments are independent ot time. 


fix - E[X(t)\ 


— E[A(t) iWuJut B{t)s\njj { )t\ 


= fl \ COS u/'()/ — flf! AnuJ()t 


For this to be. independent of t. /t A = /i B =<)=*• /i\* = 0. Similarly /* c = l'r> = 0 => fiy = 0. This 


shows that zero mean processes are required by stat.ionnrity. 


;u 


The second order moments must also he independent of time. 

Rx{t) 

- E[X(t)X(t-r)] 

=r. E[(A{t) cos it-'o t B(t) Killed) 

(A(t 7)coslJ()(/ t) - D(t r)sh\u;u(f r))] 

= R.\(t) (*os uJ\)t cos - t) + .sin o;o/ sin — t) 

R .\b(t) c< > s j t sii 1 -l’o ( t r) /?,ud r) siiiutyf awu/uf/ t) 

( 38 ) 

— (r) [cos a;nr 4 eosx*o(2t r)\ 

+^J?£?(t)[c««u;ot - cosu^(2t r)j 
^/? / \z?(" r )[ s i IliL ’ r o7‘ f siria;o(2t ")] 

— ^ /? .\z?( — ^)[ — suitor d' sin — t)] 

= f ^M^)^wW-- rl sina;or 

MlLjSalli a^,(2t r) R '** ul T) si.iu;„(2t 7 ) 

To remove the t dependence requires that 

7l.i(7-) = 7 ?c(r) and R,w(r) ~ R.\b{ t) 

{•■«)) 

=> T?.,y( 7 ) = 7 ? i(r) cosiu’ot f 77.1/3(7) sin a^ot 
A similar argument can he applied to V ( t ) to give 

7?r(r) - E d (t) and 7? t 0 ( 7 ) = 7?r/>( r) 

(10) 

-=> 7?y(r) = 7?r(r) cosj/nt f T?r/>(7) sin w it t 



The (.toss- term between the components must also be independent, of t. 


Rxy(t) 

- E[X{t)Y{t-r)} 

— E{(A(t)('osjJ{)t B(t) xmu.\)t) 

( C(t 7") COS t) D(t T ) siliwc’oft 7-))] 

— Rac{t) cos u/'o t cos (t - T ) + R r>D ( T ) s hi u,’n t sin (/■ ~ r ) 

R.j\D ( 7”) COS sill ix 1 ^ (t 7") RftC * ( 7”) sill lC{)t ('OSu/o(f 7"J 

— ^i?. 4 c( r )[cosu; u T I (tosafn(2f r)] 

f ^/?^p(r)[cos o-’oT cosu. ? (}(2f — r)] 

5^io( r )[ s hiu.\ir f sinct- ? o(2A 7")] 

- t;Rbc'(t ) [ — sin uJ\)T T sinu. i n(2f — t)] 

^ RAciry^n^r). C()S ^ jT + siu W()T 

co*uJ (] (2t r) (r - ! sin^o(2/ t) 

To remove the f dependence requires that 

Rac(t) — Rbd(t) and Rad{t) “ %’( r ) 

=> i?xv(r) - Rac( t ) cos^u/ + -R,id(t) sinu;nt 
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Appendix B Proof that the Spectral Matrix S is positive 


semi-definite 


In order to show that the Hmuitian matrix 8 is positive semi-definite it is necessary to prove the 


following i nequ ; di tv. 


< Sx.x >= x*Sx > 0 


( 43 ) 


Let ;r — [». b.c. rf). 


x'Sx 


Sa 

Sab 

Sac 

Sad 

Sad 

So 

Sn< 1 

Sbd 

S’ac 

She 

Sc 

Sod 

Sad 

StiD 

S ’ 4 

°( 'D 

So 


( 44 ) 
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Writing the S matrix in terms of Sp{uj), Sq(^j) and Spq(u>) (using Equation 2d and 2(i). 
:v*7Sx - \{Sp(u,')[oa* + hh* + i(a*b - ah*)] 


+ Sp{- 


+ w/ 1 

i(<tb* 

,rh)} 

■f-Sy -f 

M* 4 /'( 

r'd-c 

*)\ 

+S Q (- 

-OK 

4 rW* 4 

Hc<r - 

<•*(01 

+Spq{ 

ju;)[«V 

I k*d+, 

/(aM 

fc*U] 


:**)K 

1 ftrf* 4 , 

/.(/*•* 

«</*)] 

\S PQ l 

' u;)[f/r 

* 4 /«/• 

f ?!((///' 


+*7'gl 

-0[»* 

r 4 b*d 

} i{b* v 




:*. r ) 


Lot t — a 4- ib. a — r 4- id. v = a- ib. w = r - id. 

.r'lSj- = ^{5 P (u?)/r 4 5 P ( jj)rr* f Sg(u;)im* 4 S Q ( u:)irw* 

+S PQ (u)t*u 4 4 S PQ (~j;)vu' m 4 

- \{Sp(ij)tt* 4 Sp( U.')vn* 4 Sq ( ur)uir 4 Sq( 

I -2R«[S>g(-;)r«] 4 2Ro[5 Pg ( ^)r<r*]} 

> ^{5/>(u;)|fp + iSX-xOhi 2 4 $ q (*')\ u \ 2 4 Ag( ^)H 2 

2|Spg(u?)||t|M 2|.S>g( - -^)!MI W, I} 

> \{S P (*)\t\* 4 4 %MM 2 4- 5g(-u;)|w1 2 

- 2 \/ s p{oJ)Sq(ju) 1 1 1 1 f / 1 - 'l\f Sp{ - u.' ) .S q ( ~u> ) ; <' 1 1 U ' | } 

“ 4 { ( \/ *S’p(^)|^| \f Sq (uJj|u|)“ 

4 ( \J Sp{~ u/)|r| "■ \Z'S( fc ;(-^)| 7/, |) 2 } 

> o 

This proves Equation 13 and honee shows that S is positive semi-definite. 
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Figure 1: Example of realized , spectra (plotted with desired spectra) for complex processes P(t) 
and Q(t): (a) power spectral density of P(t): (b) power spectral density of Q(t): (c) complex cross 
spectra of P(t) and Q(t) (real part); (d) complex cross spectra of P(t) and Q(t) (imaginary part). 
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Figure 1: Example of realized spectra (plotted with desired spectra) for complex processes P(t) and Q(t): 
(a) power spectral density of P(t); (b) power spectral density of Q(t); (c) complex cross spectra of P(t) 
and Q(t) (real part); (d) complex cross spectra of P(t) and Q(t) (imaginary part) 











